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We calculate the current-voltage characteristic of metallic nanotubes lying on a substrate at high 
bias voltage showing that a bottleneck exists for short nanotubes in contrast to large ones. We at- 
tribute this to a redistribution of lower-lying acoustic phonons caused by phonon-phonon scattering 
with hot optical phonons. The current-voltage characteristic and the electron and phonon distri- 
bution functions are derived analytically, and serve to obtain in a self-contained way the frequency 
shift and line broadening of the zone-center optical phonons due to the electron-phonon coupling at 
high bias. We obtain a positive offset on the zero bias shift and no broadening of the optical phonon 
mode at very high voltages, in agreement with recent experiments. 

PACS numbers: 63.22.Gh, 78.30.Jw, 73.63.Fg, 73.50.Fq 



I. INTRODUCTION 

Carbon nanotubes are one of the strongest and stiffcst 
materials which can sustain very high currents before 
breaking. This electric property makes metallic nan- 
otubes an interesting alternative to nanometer-sized 
metallic wires. Since nanotubes can behave like semicon- 
ductors, their possible use in logic electronic circuits are 
promising. This has recently led to a number of exper- 
iments evaluating their current vs voltage characteristic 
at high bias voltage [H-Q, with related theoretical work 
in Ref . EME 

At low voltage, the current-voltage characteristic is 
mainly influenced by acoustic phonons and by impurity 
scattering. At higher voltage, optical phonons become 
important. For metallic nanotubes on a substrate, the 
current vs voltage curve is increasing, in contrast to sus- 
pended nanotubes where the characteristic shows a neg- 
ative differential conductivity at high bias 

We shall review in Sect. II the current-voltage charac- 
teristic of metallic nanotubes lying on a substrate. Fol- 
lowing Refs. [ll|, [lU we use a Boltzmann approach for 
the electrons coupled to zone-center and zone-boundary 
optical phonons. We take into account explicitly the dy- 
namics of the phonons by a Boltzmann equation con- 
taining an inelastic term to describe the decay of optical 
phonons into underlying acoustic phonons [ll|, [12 • We 
use first the so-called single-mode relaxation time ap- 
proximation for the scattering term jlfjj . This is char- 
acterized by a thermal phonon relaxation time r op . For 
the electron-phonon relaxation time r ep we use a numer- 
ically determined value [HI E3: which reproduce very 
well the experimentally determined lifetimes of optical 
phonons. This proceeding agrees with the numerical 
work of Refs. [ll|, [H m which the current- voltage char- 
acteristic of short nanotubes with lengths smaller than 
1 jjxa is calculated. In contrast to this, Sundqvist el 



al. (l3j have in their calculation r cp -values which are 
around three times smaller than the experimental val- 
ues. They determine the current- voltage characteristic of 
nanotubes larger than 1 /an using the one- valley approx- 
imation for the electrons. This implies that the electrons 
are scattered only by one type of phonons, i.e. zone- 
center phonons, between the bands within this valley. 
By using the experimentally determined thermal phonon 
relaxation life-time of r op « l.l±0.2ps [U|,[lj|, we reach 
a good agreement with the experimentally determined 
current-voltage characteristics of large nanotubes. This 
is in contrast with what happens in short nanotubes, 
which one has to use at least a five-times larger thermal 
phonon relaxation time to find a reasonable agreement 
with experiment. 

Due to the low dimensionality of a carbon nanotubc 
system in which a fast initial decay of optical phonons 
is followed by a slow decay of only a small amount of 
secondary acoustic phonons [20j . we expect a bottle- 
neck in the relaxation path for the hot optical phonons 
generated b y c harge carrier scattering. This idea was 
used in Ref. [2fj to explain the large discrepancy between 
the radial breathing mode lifetimes measured by Raman 
scattering experiments, and by tunneling experiments. 
Such a bottleneck leads of course to larger effective ther- 
mal relaxation times for the optical phonons. In Sec- 
tion III we shall describe this fact effectively by taking 
into account in the phonon Boltzmann description the 
secondary acoustic phonons in a simple model. By us- 
ing suitable secondary phonon relaxation times, we were 
able to reproduce the experimentally determined current- 
voltage curve also for short tubes. Our result shows that 
for long tubes the system does not exhibit a phonon 
bottleneck, in contrast to short nanotubes. We explain 
this by the fact that at tube length smaller than 1/im, 
the thermal scattering lengths of many acoustic phonons 
reaches the systems size which then dynamically closes 
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relaxation paths for the optical phonons. 

We find a similar effect in the interaction of phonons 
with the electron system under bias voltage. It was 
shown in Refs. [ll|, for short nanotubes that one finds 
a large increase of the phonon distribution function at 
the boundaries of the tube. We observe an even worse 
situation, that we do not find any numerical solution for 
the Boltzmann equation when setting the optical phonon 
velocities to zero. In contrast to this, we see for long 
tubes a phonon distribution function which is peaked in 
the center of the nanotube, in agreement with experi- 
ments for large suspended tubes [21[ . In order to under- 
stand this effect better, we solve in Section IV the system 
of Boltzmann equations for the charge carriers and the 
phonons analytically within certain approximations. We 
succeeded in reproducing especially well the large- voltage 
small-length regime of the numerical determined current- 
voltage curves. Our calculation shows that the reason for 
the increase of the temperature at the boundary of tube is 
again based on the fact that using the phonon relaxation 
path for small tubes, the electron phonon coupling part 
in the phonon Boltzmann equation creates effectively an 
additional phonon relaxation term with a negative sign. 
This leads to the increase of the phonon temperatures at 
the boundaries of small tubes. 

In the analytical calculations of Section IV we deter- 
mine the electron distribution function under high bias 
voltages. The appearance of this distribution function is 
of course much different from the Fermi distribution func- 
tion in thermodynamical equilibrium. The knowledge of 
this function opens up a number of possible applications. 
For example, in Section V, we calculate the level broaden- 
ing and frequency shift of the zone-center optical phonons 
mediated by the electron-phonon interaction under high 
bias voltage. We find a positive frequency offset on the 
zero bias shift for very large bias voltages, in contrast 
to the frequency shift mediated by the phonon-phonon 
scattering with phonons in thermal equilibrium. For very 
large nanotubes we obtain also a negative frequency off- 
set due to the electron-phonon interaction. The electron- 
phonon mediated zero-bias broadening of the zone-center 
optical mode vanishes at high voltages. The results for 
very high voltages are in agreement with a recent exper- 
iment measuring the influence of the high bias on the 
honon modes of carbon nanotubes lying on a substrate 



Summarizing, in Sect. II we discuss the results of the 
coupled electron-phonon Boltzmann system in the relax- 
ation time approximation numerically. In Section III the 
secondary acoustic phonons in the Boltzmann equation 
ia taken into account. We shall carry out in Section IV an 
analytical calculation of the current-voltage characteris- 
tic, the electron and the phonon distribution functions. 
These functions are used in Section V to calculate the 
level broadening and frequency shift of the optical zone- 
center phonons mediated by the electron-phonon inter- 
action under high bias voltage. 



II. CURRENT- VOLTAGE CHARACTERISTICS 
OF CARBON NANOTUBES 

The method we use here to calculate the current- 
voltage characteristic of metallic nanotubes is based 
on the semi-classical Boltzmann equation. Within this 
method quantum interference corrections to the conduc- 
tivity are not taken into account [22j|. It was shown just 
recently through numerical calculations that these cor- 
rections to the conductivity are negligible above room 
temperature for single-walled carbon nanotubes without 
structural defects due to phonon scattering decoherence 
mechanisms 23[. Other works using the semi-classical 
Boltzmann equation for electron or phonon transport not 
mentioned yet are found in Refs. H3, US 

The energy levels of electrons in a nanotube consists of 
one-dimensional bands positioned in the graphene Bril- 
louin zone around the K and K' points. For metallic 
nanotubes two energy bands corresponding to right (R) 
and left (L) moving electrons cross at these points. In 
the following, we assume here that the diameter D of 
the nanotube and the applied bias voltage U is so small 
that we can neglect electron excitations to higher bands. 
For example, this is valid for a nanotube with diameter 
D w 2nm when we apply a bias voltage of less similar to 
U < 2V. The electron distribution functions for a nan- 
otube under bias voltage around the K and K' points are 
equal which we denote by fL/R.(k,x,t). At larger volt- 
ages only the optical phonons are relevant as a source 
of electron-hopping between the bands. The hopping 
between one band at K and the other at the K'-point 
are mediated by zone-boundary optical phonons where 
only Kekule type of lattice distortions couple to the elec- 
tronic system [2g|. We denote in the following the cor- 
responding phonon distribution function by n K (k, x,t). 
The hopping between bands at the same K or K' points 
are mediated by longitudinal zone-center optical phonons 
with phonon distribution n r (k, x,t) [Toj . The time evo- 
lution of the electrons are governed by the semi-classical 
Boltzmann-cquation 
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d t T v F d x + —d k f L/R = [d t f L /R 



(1) 



where the collision term [dt Il/r] c ~ [dtfh/R]e + 
+ [^t/i/fl]fs + [dtfL/R\bs consists of an clastic scatter- 
ing term [d t f L ] e = v F /l e [f R (k) - f L {k)] due to acoustic 
phonon scattering (in the quasi-elastic limit) and impu- 
rity scattering. l e is the elastic scattering mean free path. 
We assume here l e = 1600 nm d, |Tl|. The electron ve- 
locity vf is given by vf = 8.4 x 10 7 cm/s and e (e > 0) 
is the electronic charge. 

[<9*/i/i?.]fs is a forward scattering term which should 
have a minor effect especially for higher applied voltages 
since it does not change the propagation direction. The 
time evolution of optical phonons is given by 



[d t + v» op (k)d x ] n v = [8 t n\ + [d t n»] 



(2) 
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where v = r, K denotes zone-center or zone-boundary 
phonons. We use here i£ p = sign(fc) 2.9 x 10 5 cm/s and 
v* = sign(fc)7.2 x 10 5 cm/s [HI, [T^ ■ The term [d t n v ] c 
is due to phonon-electron scattering, while the term 
\dtn u ] QSC represents thermal phonon relaxation. Note 
that the coupled electron-phonon system is not heated 
up by applying large voltages on the nanotube due to 
this term, which accounts effectively for the scattering of 
optical phonons into underlying (acoustical) phonons. 

Scattering of phonons with electrons leads to two scat- 
tering contributions in the electronic Boltzmann equation 
([1]) as well as in the phononic Boltzmann equation ([2]). 
When restricting on the backward scattering contribu- 
tions we obtain for the electronic scattering term 
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[n"(k + ,x) + l]f R (k R (e+))[l - f L (k L (e))} 

-n v (k + ,x)[l-f R (k R (e + ))]f L (k L (e)) 
+ n»(-k-,x)f R (k R (e-))[l-f L (k L (e))] 

- [1 + n v (-k-,x)][l - f R (k R (e-))]f L (k L (e)) 



with k^ = k R {e^) — fc/,(e) and e ± = e ± huj v . The corre- 
sponding phononic scattering term results in 



x {f R (k+)[l - f L {kl)] + h{-k- L )[l - f R (-k+)}} 

— n v {k, x)x 

x {f L {k- L )[\ - /*(*+)] + f R (-k+)[l f L {-ki)]} 
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where k^ L = k R / L {e{k/2) ± hxi u /2). The number s u is 

given by s K = 1 for zone-boundary phonons and s r = 2 
for zone-center ones (llj . In order to derive these num- 
bers one has to take into account that momentum phase 
space of the phonons is twice as large as the phase space 
of the electrons. Further one has to consider the fact 
that the electron jumps from the K to the K'-band are 
mediated by K-phonons but the reverse jumps by K'- 
phonons. On the other hand jumps of electrons within 
the same valley are mediated by the same r-phonons. 
Finally, we mention here that we used the boundary con- 
ditions FlOi 



f L (k L (e),L) 
n u (k > 0,0) 



fn(k R ((.),0) 
n v {k < 0,L) 



n F (e) . 

op 
n B ■ 



(5) 
(6) 



where np(e) is the Fermi function for a metallic nanotube 
with zero gate voltage at room temperature, i.e. np(e) — 
1/(1 + e e / fefjT ), and is the Boltzmann factor for opti- 
cal phonons at room temperature given by n°g w 0.0014. 



We use optical phonon frequencies hu K = 161 meV and 
hui r = 196 meV. The electron-phonon scattering times 
for zone-center and zone-boundary optical phonons are 
given by r r p = 538 fs and = 219 fs [l(| where we as- 
sume tube diameters of around 2.0nm typical in existing 
current-voltage experiments in the literature. 

Our method to solve ([1]) and @ is based on the numer- 
ical time integration by the standard splitting method 
[27| . We discretize the differential equations in momen- 
tum and position space [28| . To integrate the collisionless 
free electron and phonon equations in some time step, we 
use the exact solution of the equations in the case of the 
electron motion. This means that the time step value are 
fixed by the space grid. The free phonon motion in one 
time-step is given by the up integration of the collision- 
less discrete version of ((2|) on the space grid. 

In this section, we use a standard single-mode relax- 
ation time approximation for the optical phonon scatter- 
ing term [dtn"] osc given by 
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Note that this approximation is only valid for the system 
lying on a substrate. For the suspended nanotube system 
one has to take into account explictly the heat transfer 
by acoustic phonons to the leads @. We assume in our 
calculation that the thermal relaxation time T op is similar 
for zone-center and zone-boundary optical phonons. This 
approximation is justified for graphene in Ref. [13 where 
it is shown that the relaxation times of both phonon types 
are almost equal for acoustic phonon temperatures a little 
higher than the room temperature. We do not expect a 
difference for carbon nanotubes. These temperatures are 
immediately reached at the high voltage experiments we 
are interested in [|, 0, 1, [H|. 

First we calculate the current- voltage characteristic for 
nanotubes of length larger than 1/rai in the vicinity of 
T op = l.lps. This value is chosen since Song et al. [l8j 
determined experimentally r op = l.l±0.2ps in agreement 
to the experiment of Kang et al. in Ref. |l9|. 

In the upper panel in Fig. 1 the current-voltage char- 
acteristic determined with help of Eqs. ([I])-® is shown 
for nanotubes at bias voltage U = EL = IV as a func- 
tion of their length for various relaxation times r op in 
the vicinity of the experimentally determined relaxation 
time. The solid curve in the Figure is given by the exper- 
iment carried out by Sundqvist et al. in Ref. 0. In this 
experiment the nanotube length was effectively varricd 
by changing the distance between the electrodes where 
the bias voltage is applied. 

We get the best agreement within the experimental 
uncertainties for the quantity r op = 1.1 ± 0.2ps at value 
Top = 0.9ps. In the lower panel in Fig. 1, we calculate 
the current- voltage characteristic of nanotubes for r op = 
0.9ps as a function of the tube length for various bias 
voltages. Fig. 2 shows the average phonon density for a 
nanotube at length L = 3000nm, U = IV and r op = 0.9. 



This density is determined by 



(x) 



2eEL 



eE(x-L) 



eEx 

de [n v (2k L (e),x) + n v (2k R (e),x)} 

(8) 

The factor two in the denominator is necessary due to 
the fact that we average over the right and left moving 
electron bands. We obtain phonon densities which are 
peaked in the center of the nanotubes. This behaviour 
is in accordance to experiments [2l{ | for suspended nan- 
otubes. We show also in this figure the space averaged 
phonon densities n v = J dxri v (x) /L as a function of the 
tube length for various bias voltages. 

One reason for the small difference in the current- 
voltage characteristic between experiment and theory in 
the upper panel of Fig. 1 at large nanotube lengths is 
due to the fact that the diameter D of the tube in the 
experiment Ref. is in fact a little larger than 2nm. 
Although the diameter was not measured explicitly in 
Ref. [H one can estimate it by the fact that Sundqvist et 
al. measured approximately half of the differential resis- 
tivity for short distances between the elctrodes in com- 
parison to the value in Refs. @, 0|- In these experiments 
the current- voltage characteristic of short nanotubes with 
a measured diameter D sa 2nm was recorded. By using 
an analytical theory for the current- voltage characteris- 
tic which will be derived in Sect. IV and further that 
the electron-phonon scattering time t" is proportional 
to the diameter D of the nanotube [1(| we obtain that 
D w \pl x 2nm w 2.8nm. We show in the upper panel 
in Fig. 1 by the dotted curve the theoretically calculated 
current-voltage characteristic for a 2.8nm nanotube, i.e. 
Tg p is now a factor v2 larger than for the 2.0nm nanotube 
already used before, and r op = 0.9ps. 

In contrast to the small undershooting of the theoreti- 
cally determined current- voltage curve in comparison to 
the experimental curve for large lenths in Fig. 1, we ob- 
tain for small nanotube lengths an overshooting of the 
curve. The reason for this different behavior between 
large and small nanotube lengths will be discussed in the 
following. 

Next, we determine the current-voltage characteristic 
of short nanotubes. In Fig. 3 we show the current-voltage 
characteristic for a nanotube of length L = 300nm as 
a function of the bias voltage U for various phonon re- 
laxation times r op . The solid curve is given by the ex- 
periment [||. We obtain the best agreement between 
experiment and theory for the differential conductivity 
dU/dl « 220fcSl at r op « 9.1ps. It is astonishing that 
this optical phonon relaxation time is much larger than 
the experimentally determined optical phonon relaxation 
time r op = 1.1 ± 0.2ps. The reason for this discrepancy 
will be discussed in the next section. Note that we ob- 
tain in Fig. 3 in the low voltage regime a better agree- 
ment between the theoretically and experimentally deter- 
mined curves by using smaller elastic scattering lengths 
L < 1600nm. 
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FIG. 1: (Color online) Upper panel shows the current- voltage 
characteristic of long nanotubes with diameter D — 2.0nm 
for bias voltage U — IV calculated by the help of H}-© 
for various thermal relaxation times r op (solid curves) and 
l e — 1600nm. The dotted curve shows the current- volt age 
characteristic of a D — 2.8nm nanotube, i.e. it uses y2Te P as 
the electron-phonon scattering times, and r op — 0.9ps. The 
(black) circles are given by the experiment [sj. The lower 
panel shows the current-voltage characteristic of long nan- 
otubes with diameter D = 2.0nm for various bias voltages, 
Top = 0.9ps and l e = 1600nm. 



III. SECOND GENERATION PHONONS 

From Fig. 3 we see that a satisfactory agreement be- 
tween the experimentally and numerically determined 
current-voltage characteristic is only reached for T op S> 
l.lps. On the other hand, recent phonon lifetime 
experiments on carbon nanotubes show that t od ~ 

■ ■ I I " 

l.lps for zone-center phonons [18j, llSj]- These mea- 
sured phonon lifetimes are governed by the decay of 
zone-center phonons to two lower energetic second- 
generation phonons where the number of these decay 
channels should be rather small for one-dimensional nan- 
otube systems in contrast to higher dimensional systems 
like graphene or graphite [20| . The second-generation 
phonons are typically acoustic ones which then again 
scatter in two acoustic phonons with even lower energy 
and longer wave-length where this lifetime is much longer 
than of the primary optical phonons. The reason for the 
longer lifetime comes from the fact that the three phonon 
matrix element vanishes in the long- wavelength limit and 
further that the phase space for phonon decay is smaller 
for lower phonon energies due to energy conservation. 
The long lifetime of secondary phonons and the small 
amount of possible decay channels could lead to a bottle- 
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FIG. 2: (Color online) Left panel shows the phonon-density 
distribution functions n K (x) (solid curves) andn r (i) (dashed 
curves) defined in ([EJ for U = IV, r op = 0.9ps and L — 
3000nm. The black curve is calculated by using the former 
defined «£, = sign(fc) 2.9 x I0 5 cm/s and v^ p = sign(fc) 7.2 x 
I0 5 cm/s [HLG3I- The red curves uses «o P /5 and v^ p /5 as 
phonon velocities. The curves lie practical on top of each 
other. The right panel shows position averaged phonon dis- 
tribution functions n K (solid curves) and rf (dashed curves) 
for U — 0.4V (lower red curves) and U — IV (upper black 
curves) for r op = 0.9ps. 



, l.lps 




FIG. 3: (Color online) Current-voltage characteristic of a 
L — 300nm nanotube calculated by the help of for 
various thermal relaxation times r op (dashed curves) and 
l e = 1600nm. The (black) solid curve is given by the ex- 
periment [J]. The (red) dotted curves are calculated for 



9.1ps, l e — 800nm (upper curve) and Z e = 400nm (lower 



curve I 



neck in the decay process. This means that a significant 
amount of secondary phonons are assembled in the decay 
of hot-phonons generated by charge carriers through the 
electron-phonon interaction. When this non-equilibrium 
amount of secondary phonons is similar to the number 
of equilibrium phonons following the Bose-Einstein dis- 
tribution the single-mode relaxation time method lead- 
ing to the appearance of the scattering expression (J7J is 
no longer valid. In Ref. it was argued that this fact 
is responsible for the considerable difference in the life- 
time measurements of the radial breathing mode by using 
either Raman-scattering experiments or electron tunnel 
experiments. The decay channel can then be described 
by the following Boltzmann equations [3ll ] when neglect- 
ing the phonon velocities v^ p w on the left-hand side 
of ©: 



d t n v = — \-n v {l + n ac ) 2 + (1 + n")nU + [d t n v 



(9) 



P^op 



K(l+n ac ) 2 -(l+„>2 



— (n-ac - ?4 C ) 

7~a.r. 



(10) 



Here the first term in the brackets in © describes the 
scattering of the optical phonons with distribution func- 
tion n v into two secondary phonons with distribution 
function n ac . For simplicity we assumed that the sec- 
ondary phonons follow all the same distribution function. 

The second term in the brackets in ([9]) describes the 
reverse process. The second equation (fTt)|) describes the 
dynamics of the secondary phonons. Here p denotes the 
number of decay channels. For one-dimensional solids 
this number is generally small [2Cj| . For simplicity we 
further assumed in (|10j) that the secondary phonons are 
coupled to a heat bath where the relaxation with this 
bath happens in time r ac 3> r op . The quantity is 
the Bose-factor for the secondary phonons which we as- 
sume to be half of the frequency of the optical primary 
phonons, leading to n|f ps 0.03 at room temperature. 
Note that this choice is consistent with the fact that we 
choose uniform secondary phonon distributions. 

In the stationary case we have dtn ac = 0. Then we 
can solve the second equation (fT0|) for n ac and insert the 
result into the first equation which leads to an effective 
optical phonon scattering term 
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FIG. 4: (Color online) Upper panel shows the current- voltage 
characteristic of a L = 300nm nanotube calculated by the help 
of CP"© an d pip for various parameters Tac/piop (dashed 
curves) and r op = l.lps, l e = 1600nm. The (black) solid 
curve is given by the experiment [j. The lower panel shows 
the current-voltage characteristic of various nanotubes with 
different lengths for r op = l.lps and r a c/pTop = 20, l e = 
1600nm. 
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FIG. 5: (Color online) Left panel shows the phonon-density 
distribution functions n K (x) (solid curves) and n r (x) (dashed 
curves) defined in |(8j for U = IV and r ac /pT op = 20, 
r p = l.lps and L — 300nm. The black curve is calcu- 
lated by using the former defined phonon velocities Vp P = 
sign(fc) 2.9x 10 5 cm/s and = sign(fc) 7.2x 10 5 cm/s [Tl|,E3. 
The red curves uses t>o P /5 and Va P /5 as phonon velocities. The 
right panel shows position averaged phonon distribution func- 
tions n r (dashed curves) and n K (solid curves) for U = IV, 
T a c/pTop = 20 and L = lOOnm (blue curves), 300nm (black 
curves) and 500nm (red curves). 
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PT op /T^ c ^0 
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As is seen from (fT2"|) in the case of no existent bottleneck, 
i.e. large pT op /T ac and small Boltzmami factors n°g, rig 
valid in our case, we obtain the standard single-mode 
relaxation time approximation ([7|) for the optical phonon 
scattering term. 

In the following, we carry out the numerical calcula- 
tion by using (JTJ)— (J3J) with the optical phonon scatter- 
ing term (fTTj) substituting (JT]). From (fTTj) we obtain 
that the optical scattering term depends via T a c/pT p 
on the acoustic scattering length. We show in the up- 
per panel of Fig. 4 the current voltage characteristic 
for various parameters T ac /pr op , r op = l.lps and nan- 
otubes of length L ?s 300nm. This figure should be 
compared to Fig. 3 in the case of the single-mode relax- 
ation time approximation which uses (J7]) for the optical 



phonon scattering term. We obtain a similar behavior of 
the current-voltage characteristic curves in both approx- 
imations. The reason is seen from expression (|13[) . Wc 
obtain from this expression and the fact that at room 
temperature n^,n°^ < 1 as well as > 1 for high 
bias voltages, that for small pT op /r ac the effective optical 
phonon scattering term [dtn v \ osc has still the standard 
single mode relaxation form ([7]). The effective relaxation 
time T op is then changed to r op — > r ac /2p which respects 
the fact that the relaxation of the optical phonons are 
effectively relaxed on 2p channels of relaxation time r ac . 

In the lower panel of Fig. 4 we show the current- voltage 
characteristic for various nanotube lengths by using (fTTj) 
as the optical phonon scattering term with T ac /pr op = 20 
and T op = 1.1 ps. Fig. 5 shows the energy-averaged 
phonon distribution function ([5]) n v (x) (left panel) and 
the energy and position-averaged phonon distribution 
function rf (right panel) for various nanotube lengths. 
We obtain that in contrast to the case of large nanotubes 
in Fig. 2, the short nanotubes show an increasin g ph onon 
density at the boundary of the nanotube [ll], EH ■ We 
obtain from the left panel in Fig. 5 the even worse fact 
of a diverging current-voltage characteristic for phonon 



velocities 



op 



0. We shall understand this unusual 
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behaviour better in the next section where we show an 
analytic solution of the Boltzmann system (JT])-([7|). 

Summarizing, by taking into account a possible bottle- 
neck in the optical phonon relaxation path we obtain for 
an optical phonon relaxation time r op = l.lps for large 



nanotubes an effective scattering parameter T ac /pr 







as the best fitting parameter to the experimental curves, 
i.e. no phonon bottleneck is seen in this case. On the 
other hand for small nanotubes and r op = l.lps we ob- 
tain an effective scattering parameter of T ac /pr op = 20. 
The reason for this discrepancy between large and small 
nanotubes lies in the fact that we have neglected the 
velocities of the secondary phonons in the Boltzmann 
equation (|10p . This leads in the case of the decay of 
zone-center optical phonons to additional terms of the 
form ±v ac d x n^ c in the left hand side of (fTU|) where nf c 
stands for the left and right moving phonon in a scatter- 
ing pair. This means that we remove the restriction in 
the right hand sides of ([9]) and (fT0|) that the scattering 
pairs have all the same distribution function. To see when 
the ivacdxTiac term becomes relevant in the Boltzmann 
equation we further have to take into account the num- 
ber of decay channels p for nanotubes. In the following 
we restrict ourselves to the relaxation of the zone-center 
phonons. The argument for the zone boundary phonons 
works similar. 

It was shown in Refs . [20L l30l that in the case of graphene 
the zone-center optical phonons scatter into three sorts of 
different pairs of phonons lying on rings in the Brioullin 
zone around the T points where scattering into the lon- 
gitudinal acoustic sector is in fact the most dominant. In 
order to estimate from this fact the number of pairs for 
a nanotube with diameter of around 2nm we use in the 
following the zone folding approximation method. By 
using as an approximation that in the case of graphene 
the decay rings lie in the mids between the T and K 
point (best fulfilled for the longitudinal mode [2(], |30| ) 
we obtain as an estimate for the number of decay pairs 
p < 20 in a 2nm nanotube. In order to see no bot- 
tleneck in the relaxation process for large nanotubes we 
have T ac < PT Q p meaning that r ac < 22ps. On the other 
hand for a nanotube of diameter 2nm one obtains for 
the lowest lying phonon modes in a simple model relax- 
ation times which are larger than around r ac > 20ps [32| . 
This leads us to the estimate r ac ~ 20ps for the effec- 
tive acoustic relaxation time. With v ac ~ 21km/s (we 
choose the maximum velocity value for acoustic phonons 
in graphite j33l - l35j ) we obtain that the acoustic phonon 
velocity terms ±v ac d x n^ c ~ ±v ac n^ c /L in the phonon 
Boltzmann equation becomes relevant for L < 420nm. 
This is only a very rough approximation for this length. 

At this length we find that for one participant of the 
scattered acoustic phonon pairs this additional relaxation 
term is not relevant since the relaxation path is already 
open. For the other participant this term leads to a clos- 
ing of the relaxation path which we saw in our numer- 
ics as an increase of the effective relaxation parameter 

7"ac / PT~op ■ 



IV. ANALYTIC CALCULATION 

Due to the similarity of the phonon frequency of zone 
boundary phonons lj k and zone-center phonons w r and 
since the electron- phonon coupling constants s k /t^ and 
s r /r^p are similar we use in the following the simpli- 
fication that the nanotube system interacts with only 
one sort of phonons with frequencies to. The effective 
electron-phonon scattering parameter r cp in the elec- 
tronic sector, the electron-phonon scattering parameter 
s p /r ep in the phononic sector and phonon velocities i; p h 
are chosen in the following way 
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sgn(k) 5.95 x 10 



(14) 



In the following discussion, we use the abbreviation = 
vft cp / (2n + 1) for the reduced effective scattering length 
and l sc = L/(l + L/ll c ) for the total scattering length. 

Below, we solve the Boltzmann equations ([I]) and ((2|) 
for large voltages ell ^s> huj and lengths L ^> ll c ana- 
lytically by the help of two approximations. In the first 
approximation we use in the electronic Boltzmann equa- 
tion (fl]) positional and momentum independent phonon 
distribution functions n(fc, x) « n. We shall determine n 
then similarly to ([8]) where the energy average is taken 
over those energies where n(fc, x) or /l(A;, x), fn(k, x) are 
non-zero, respectively. Thus we neglect large (infinite) 
energy regions where the electron and phonon distribu- 
tion functions are zero since they do not contribute to 
the current. 

In the calculation below, we find for the energy aver- 
aged phonon distribution functions 



n(x) 



where 



and 



1 



2(e u - e d ) 



de n(2kL(e), x) + n(2fc^(e), x) 

(15) 



eEx for n 
for 7i 

7T/ SC 



eUl B 
.ebl 



> flUJ, 



eE(x - L) for > hu , 







for 



.eUl s . 



n is then determined by the average 



If 

j dx n(x) . 
L Jo 



(16) 



(17) 



(18) 



The second approximation is given by a linearization of 
the non-linear scattering terms in ([3]) which can be iden- 
tified by extracting the brackets in (J3J). These terms are 
equal to the terms followed by setting in © n" equal to 
zero. To linearize these terms we should take care on the 
expansion points f L where fi, = f L + A/l and similar 
for f R . In a first crude approximation we use in the fol- 
lowing as the expansion point f L = f R = 1 which are the 
boundary values ([5]) for fi, and f R on the non-zero mo- 
mentum support of the electronic distribution function. 
Then we obtain 

f R (k R (e + Huj v ))[l-f L (k L (e))} 
- fL(k L (e)){l - f R (k R (e - fuv»))} 

~ f R (k R (e-hw»))- f L {k L (e)). (19) 
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FIG. 6: (Color online) We show the function Bi(x) defined 
in (TA13| . 



A. Current-voltage characteristic and electron 
distribution function 

With the help of the approximations (fil^) and (flU]) one 
can solve §T§ with © using Fourier methods. After a 
lengthy calculation carried out in Appendix A we ob- 
tain for the electron distribution function /l(x, kr 1 (e)) 
(|X5j) with (fX6|) and (jAlljl . The distribution function 
f R (x,k R (e)) is then given by /l(x, fci(e)) with the help 
of the substitution (|A4[) . The current voltage character- 
istic is given by (|A12[) 



HujBi 



i + L/r ac 



eU 



(20) 



where Bi(x) is defined in (|A13[) and plotted in Fig. 6. 
By taking into account the regime fi > 1, vft op ~ 130 
we obtain for L = 300nm that B\ w 0.5. By using 
l sc ~ 10.4nm (this will be shown below by using Fig. 7) 
we obtain excellent agreement with the numerically de- 
termined current-voltage characteristic at high voltages 
shown in Fig. 3 and with experiments 0, H| measuring 
l sc si 10- limn. Note that the first term in ([20]) cor- 
responds to the y-axis value obtained by extending the 
high- voltage curves to this axis. For this we further note 
that 4(e/h)hoj m 26.6/iA. 

Next, we discuss the current- voltage characteristic for 
large lengths, i.e. l sc /L <C t cp /t op and low voltages 
ireU <C HujL/I sc . Here, we obtain n w 1.2 (see the dis- 
cussion in Subscct. B.2) leading with Fig. 6 to B\ sa 1.3 
for nanotube lengths L 3000nm. This leads to cur- 
rents which are approximately two times the current val- 
ues shown in Fig. 1. Although this value is too large, 
the overall behavior of an approximatively length inde- 
pendent current for fixed bias voltage shown in the nu- 
merical calculation in Fig. 1, is also seen in the analytic 
calculation. 

Summarizing, from (|20[) we obtain a different be- 
haviour of the current voltage characteristic in the regime 
neU ^> fnuL/l sc where the second term in (|2U1) is the lead- 
ing contribution to the current and the regime irelJ <C 



HujL/I sc where the first term is most relevant. The rea- 
son for this can be seen in ([1} and ([3]). At low optical 
phonon scattering rates, i.e. irelJ 3> hioL/l sc , electron 
scattering takes primarily place from the upper part of 
the filled right moving band to the empty part of left mov- 
ing band. For higher scattering rates irelJ <C fvjjL/l sc a 
large amount of electrons are also able to be scattered 
from the upper part of the filled right moving band to 
the filled part of the left moving band within many scat- 
tering processes where now Pauli blocking prohibits this 
scattering. This Pauli blocking is only roughly described 
by the linearized phonon scattering approximation (|19| . 
This is the reason that in the large length regime we ob- 
tain less agreement between our numerical and analytical 
results in contrast to nanotubes of smaller lengths. 



B. Phonon distribution functions 

Next, we determine the phonon distribution function 
n(k,x) by solving ((2]) with ((U and ([7|). We use here 
in our analytical calculation for simplicity the standard 
relaxation time approximation ([7]) for the optical scat- 
tering term instead of the more complicated scattering 
term (fTTjl which takes into account also the second gen- 
eration phonons. As was discussed in Sect. Ill the differ- 
ences in the current-voltage characteristic are only minor 
when taking the effective optical phonon relaxation times 



' op 



0.9ps for long nanotubes and r op «a 9.1ps for short 
ones. We shall determine first the phonon distribution 
function n(k,x) in the regime ireU ^> hu>L/l sc . 



Phonon distribution function for ireU 3> hujL/l s . 



We obtain in App. A for n(k,x) (|A18j) . (|A19|) and 
(|A20[) for k > in the different parameter regimes (|A15[) . 
For k < n(k, x) is given by n(—k, —x) (IA14|) where the 
Kis are defined in (|A17[) . From these equations we can 
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FIG. 7: (Color online) The (black) solid curve is given by the 
phonon distribution function n as a function of the nanotube 
length L calculated with the help of (TA"26j) with (TX27) , (1X28)) 
and (|A29|I . The (red) dashed curve shows K\ (jA21j) as a 
function the nanotube length L in the parameter regime IZ2 
(|A17|) using n. 



calculate the energy averaged phonon distribution func- 
tion n(x). This function is given by (|A22p with (|A23|1 . 
(|A24p and (j A25|) . From this function we obtain for the 
position averaged phonon distribution function n (IA26I) 
with (fA"27j) , (|A"28|) and (|A29j) . This function determines 
effectively the averaged phonon distribution function by 
using the definition for Ki (jA21j) in the regime 72-2 (|A15jl . 
We solved this equation analytically below Eq. (|A30[) in 
various nanotube length regimes. In Fig. 7 we show the 
numerical solutions for n and K\ in the regime H.2 as a 
function of the nanotube length L. For a nanotube of 
length L = 300nm we obtain / sc w 10.4 and K\ w 0. 
For these values we obtain (|A31|) for n(x) showing in 
fact a phonon distribution function which increases at 
the boundary of the nanotube. From (|A25[) we obtain 
that this behavior is even more pronounced for A'i-valucs 
larger than zero. 



2. Phonon distribution function for -k ell <C hwL/l sc 

In the following, we restrict ourselves to the regime 
where l sc /L <C t op /t op which is good fulfilled in the large 
length regime considered in Sect. II. By taking into ac- 
count (|A6j) . (jAllj) and (jA4|) we obtain that n(k, x) ^ nff 
only in the regime \it(hvF\k\/2uj)L/l sc \ < 1 where we take 
into account that n < 1 as will be shown immediately be- 
low. 

This leads to (|A34p for n(k,x) and ()A35[) for n(x), n 
in this regime. By taking into account s p t op /t cp w 4.76 
we obtain n « 1.2. Note this value is much larger than 
the numerical determined values shown in Fig. 2. 



V. THE ELECTRON-PHONON COUPLING 
INDUCED FREQUENCY SHIFT AND 
LIFETIMES OF OPTICAL PHONONS AT HIGH 
BIAS 

By using the electron distribution functions /j, (|A5p 
with (|A"6j) . (jAllj) and f R (|A"4|) we are now able to cal- 
culate the effect of a large bias voltage on the frequency 
shift and lifetimes of optical phonons. We restrict our 
calculation to the Raman active zone center T phonons. 
In this section, we shall carry out a similar calculation as 
was done in Ref. for the frequency shift and lifetimes 
of optical phonons at zero bias voltage and temperature 
T = 0. There, the electron distribution functions Jl and 
fa are Fermi-functions. The retarded phonon self-energy 
at zero momentum is given by 

1/2 

n L ' T {u) = Y^ f dt[ii L ' T (n,u)-ii L > T (n+t,o)] (21) 

" -1/2 

where H L ' T \n,ui) is the self-energy contribution of the 
?ith electron band. Its value is given by [23] 

n-(„,.> = -4££*(£) 2 ^ (22) 

s,s' k v 7 

x i( 1± s*VM-fc 2 n f[e s (n,k)}-f[e s '(n,k)} 
2 V K 2 [n] + fc 2 ) - e s (n,k) + e s '{n 1 k) + iO 

where the electron bands for metallic nanotubes are given 
by 

e s (n,k) = shv F ^Jn{n) 2 + k 2 (23) 

with n(n) = 2im/L . s = +1 for the energy levels in 
the conduction band e > and s = — 1 for the energy 
levels in the valence band e < 0. is given by the 
number of unit cells and M is the mass of a Carbon 
atom. The subtraction of the last term in ([21j) is due 
to the fact that in order to calculate the frequency shift 
for nanotubes, we insert in the calculated expressions the 
known optical frequencies of graphene which then results 
in a double countin g w hen we only use the first term in 
(j2Tjl as self-energy [23]. The valley degeneracy is here 
considered by a factor two in correspondence to similar 
expressions in Ref. H^. The upper sign corresponds to the 
self-energy of longitudinal phonons H L , the lower sign to 
the transversal ones n T 

In the case of the frequency shift of the longitudinal 
and transversal optical T mode Ishikawa et al. use the 
zero temperature Fermi- function for f[e s (n, k)]. The fre- 
quency shift Aw and broadening T is given by [29| 

Alu = Rc[n(u; r )] , T = -Im[n(u; r )] . (24) 

In Ref. [23" it is then shown that this leads to a good 
agreement of the theoretically calculated frequency shifts 
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and broadenings by using (f2T]) - (f24|) and the experimen- 
tally determined ones using Raman spectroscopic meth- 
ods. In the following, we carry out a similar calculation 
for the case of the electron system under high bias volt- 
age. For this we use for f[e s (n, k)] in ([22)) for the lowest 
energy band, i.e. n = 0, the distribution functions /l 
and fn calculated in the last section. For the higher 
bands n^Owc shall use the Fermi-function since we did 
not take into account higher band excitations in Sect. IV 
being negligible in the considered voltage regime. 



f[e s (n,k)} 



h(k) for 

/fl(fc) for 
np(e s (n,k)) for 



?i = & k-s<0 
n = & k-s>0 . 

(25) 

Expression ([22]) for the phonon self-energy corresponds 
to a current-current Green's function loop being the 
lowest order approximation for the phonon self-energy. 
When using dressed Green's functions by taking into ac- 
count the electron-phonon interaction and also the exter- 
nal electric field one has to use non-equilibrium Green's 
function techni ques in order to get the corresponding 
loop expression [36[ for the retarded phonon self-energy. 
The relevant Green's functions in the loop are given by 
lesser G < and greater Green's functions G > . We now 
express these Green's functions by the spectral function 
and use the quasi-particle approximation (36j established 
for deriving the Boltzmann equation. This leads to (f2~Tj) , 
(ggj) with (ggp. 

It is well known that, by using the free Green's func- 
tion in loops as was done by Ishikawa et al. [29 |. this 
approximation is conserved [13] which means that the 
charge-current response functions corresponding to the 
loop fulfill the continuity equation 3JJ ■ By using a more 
general dressed Green's function in the loop one needs 
also vertex corrections in order to fulfill the continuity 
equation. It is straight forward to show that the current- 
density correlation functions corresponding to ([21]) , ([22]) 
with ([25]) which consist of dressed Green's function in 
loops with the additional quasi-particle approximation 
does indeed fulfill the continuity equation. This justi- 
fies approximation (f2Tj) . (|22|) with ([25]) for the phonon 
self-energy at zero momentum to calculate the phonon 
self-energy under high bias voltage. 

In the following, we use the abbreviations: 



ui- 



a(D) 



■kD 
2itvf ' 
27 a2 ( 2(hv F ) 



7T \ 



2Ma i 



1 



a 
ttD 



(26) 
(27) 



Here a is -\/3 times the equilibrium bond length. With 
the knowledge of the parameter T^ p = 538fs determined 
by density functional methods for a nanotube of diame- 
ter 2nm [ll| we are able to determine /3 and thus a(D). 
By using that r^ p is more generally proportional to the 
diameter of the nanotube [1(| we obtain from ([2"T]) that 



a(D) w Q.la/nD. 



(28) 



With the help of (|2Tj), ([22]) and ([251) ty using the ab- 

&=o(w) + n& >0 ( 

2 



breviation H L (u>) = fl|j =0 (w) + n^ >0 (w) we obtain 



n T (w) = a(D)Lu F [ 2 - ^-w 2 
9 



(29) 



n£ =0 (w) = a(D)ui r ^]n^ - ^ 2 -i^ . (30) 

II^ >0 is a correction factor to the zero bias self-energy 
n^ =0 ([30]). By taking into account that H T (n = 0, ui) = 
([22jl we obtain that this factor is of similar order as 
the zero-bias self-energy only in the longitudinal sector 
for ell <§; hvp2/D. Note that we restrict here ourselves 
also to the regime ell <§C Hvf2/D as was done in the last 
sections by taking into account that excitations to higher 
bands are negligible. We obtain for the self-energy n^ >0 



n 



u>o 



(w,x) = (uj,x) + (w,x) 



(31) 



with 



Uf(uj, x) = a(D)Lo r \ 2 C(x) + 2 In 



2eU 
fvjj 



U^(u J ,x)^-a(D)u; r -{(l-S ; )\n' 



(32) 



1-x 




X 



-i \ 1 -2-Wzfc\ \[ , f to L hw_ 
-xm — > + < — > — n 



eUir l sc eU J ' 
(33) 



where 



C{x) = x ln(l -x) + (l-x) \n(x) . 



(34) 



and x = x/L. The last line in (|33p means that we have 
to add the foregoing expressions with the substitution 
(lj/it)L/I sc — > -wn. Here is the self-energy part cal- 
culated b y the help of f t=0 ([Ml in ([22]). 11% is the result 
for (lATTI) . In ([331 we carry out for the logarithmic 
term in the electron distribution function (|A9[) the ap- 
proximation that we set this term constant over the range 
7r |e — eExi\l sc /hujL < 1 with its value at |e — eEx±\ = 
and zero elsewhere. This approximation leads to the log- 
arithmic singularities in (|33p . They are softened when 
using the exact functions (|A11[) without approxi- 

mation but this treatment has the disadvantage that we 
would not obtain analytical results. We should addition- 
ally mention that we neglect those imaginary terms in 
([32]) and ([33]) which exist only in small regions in posi- 
tion space of length Ax = Lhuj/eU. These correspond to 
regions where the ln-terms in (|32j) and ([33]) gets singular. 

By taking into account the results in ([317]) . ([3"T]) . ([52"]) 
and ([33]) we obtain the surprising result that the imagi- 
nary part of the longitudinal phonon self-energy and thus 
the level broadening vanishes 



Im[n L (u;,x)] = 



(35) 
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in the high voltage-bias regime eU ^> hu> in contrast to 
the case of no bias (|30|) . The reason for this vanishing is 
seen in (|2"2")h We only obtain an imaginary part for the 
phonon self-energy when there is a substantial changing 
of the electron distribution function in the energy range 
\2e\ < hw. In the high bias regime the electron distribu- 
tion function (|25p becomes constant in this range which 
is not the case for the Fermi-function in the zero-bias 
system. 



A. Position averaged self-energy 



Finally we calculate the position averaged self-energy 
nf = / dxllf(u},x)/L. We obtain 



n : (w) = a(D)u l 



II 2 (w) w a{D)u l 



-3 + 2 In 



2eU 



hw L 
ireU l sc 



, hui 



(36) 
(37) 



with 

Fi(x) 



1 

4 
+(a; 



(x + l) 2 [ln(|x|) - ln(|l + x\)] - - ln(|x|) 



-x) 



(38) 



We show in Fig. 8 the function Fi(x). In the 
regime of very high voltages eUl sc /L ^> haj/n 



we obtain that il 2 (uj) is negligible in compari- 



B. Line-shape of Raman signal 

From the considerations above it is not clear how the 
actual line-shape of the Stokes or anti-Stokes signal looks 
like in an actual Raman scattering experiment. From 
above we obtain first that the Raman signal correspond- 
ing to the response on transversal phonon mode excita- 
tions denoted by G + is not changed from the zero bias 
result. This is not true for the G~ Raman mode corre- 
sponding to the response on longitudinal optical phonons. 
In order to get a better insight into the actual Raman 
line-shape we assume that the incident laser light illumi- 
nates the nanotube continuously over the whole width. In 
typical Raman experiments the scattering of the electron 



system with the incident light is dominated by the res- 
onant scattering of a valence band and conduction band 
of fixed index n > (|2"3"|) (38j . A further enhancement of 
the signal is reached when electrons from the band edges 
are scattered. This is due to the van Hove singularities of 
the density of states at this region leading not until then 
to the opportunity of measuring phonon and electronic 
properties of single nanotubes. Further let us assume in 
a first approximation that the phonon distribution func- 
tion n T (x) is homogeneous over the nanotube width 
which was also assumed in our analytical calculations in 
Sect. III. This leads us to the conclusion by taking into 
account that only the electron distribution function of 
the lowest electronic energy band n = is changed due 
to the large bias voltage that we can determine at least 
approximately the line-shape of an actual Raman signal 
by the position average over the individual Raman sig- 
nals. 

In the following, we use the abbreviation 



F 2 (x) 



dy S(C(y) - x) 



(39) 



;Fo 



son to nf(oj). For nanotubes of diameter 2nm 
we have Qq ~ 0.35. Thus we obtain that 
Re[n^ =0 (w r )] w -0A6a(D)uj r leading to the re- 
sult that Re[n L (w r )] w a(D)uj r [-3A8 + 2 ln(2e[//fiw r )]. I cp (cj)oc{^F 
This means that we obtain in the high voltage regime 
for eU/huj r > 1.8 a positive frequency shift for the 
longitudinal optical phonon frequency at the T point. 
On the other hand, by taking into account that 
lim^^oo F\(x) = — l/21n(a;) we obtain that for fixed 
eU/hu> r and in the large length limit L/l sc — > oo where 
n < 1, that the frequency shift is negative. 



in order to calculate the line-shape of the Raman signal. 
We restrict ourselves in the following to the most impor- 
tant high bias regimes ireUl sc /L 3> hu r and to the large 
length regime where ireUl sc /L <C hui r , but eU 3> hum 
where we obtain simple expression for the Raman signal. 
By taking into account that Im[II L ] = we obtain for the 
line-shape of Stokes and anti-Stokes signal / e p by using 
(EH), ([331 and (OH 



Re[n£ =0 (^ r )] ] n (2eU) 



2a(D)uj 



= o(" r )] 



for^»-^ 

huj 1 7tI sc 



(5/2)a(D)uj r 



ln( 



4 i„ I 2eU\ 

hu> T ) 



1 -KeU l t 



for 7I « 1 gL«^. 



(40) 

We show on the right hand side in Fig. 8 the function 
F-2. Due to the the maximum of C(x) at x = 1/2 with 
C(l/2) = — ln(2) we obtain that F 2 (x) is singular at 
x = — In 2. This singularity leads to a sharp edge of the 
Raman spectrum at the corresponding frequency accord- 
ing to (|40"1). 

By taking into account also the broadening of the 
phonon modes due to the phonon-phonon scattering we 
obtain that the signal of the actual Raman mode is a 
convolution of I cp (u>) (|40[) with a Lorentzian which has 
a negative frequency shift and broadening correspond- 
ing to the phonon-phonon interaction contribution to the 
phonon self-energy. 



C. Discussion 

Next we compare our findings with recent experiments. 
In Refs. IH, Hi] the optical phonon lifetimes and frequency 
shifts are measured for suspended metallic carbon nan- 
otubes at high bias voltage using Raman spectroscopy. 



These investigations did not find a response of the G_ 
and G+ modes corresponding to the longitudinal and 
transversal optical T modes for all measured nanotubes 
where only the former mode couples to the electron sys- 
tem at high bias. In fact, it is not well understood which 
mode response to the bias for a certain nanotube. The 
nanotube phonon temperature was then measured by two 
different methods either by the Stokcs/anti-Stokes in- 
tensity ratio being a function of the phonon Boltzmann 
factor of the corresponding mode or by mode softening 
where one compares the softened mode with the known 
temperature softening which is caused by phonon-phonon 
scattering. Astonishing, that both methods find phonon 
temperatures which are in good accordance. They did 
not provide any hint for a different behavior which we 
found theoretically above. 

Such a different behavior was in fact seen by Oron- 
Carl et al. @ for metallic nanotubes on a SiC>2 sub- 
strate. They found a mode hardening and a reduction of 
the linewidth at high bias voltage in comparison to the 
zero bias voltage values. A real quantitative comparison 
of our results with the experimental finding is difficult 
here since the overall broadening and hardening is the 
additive effect of phonon-phonon scattering leading to 
softening and broadening of the mode and the harden- 
ing effect and a linewidth reduction due to the electron- 
phonon coupling. That this non-thermal effect was not 
seen in free-standing nanotubes could be explained by the 
fact that the heating of many of the crystal modes are 
much easier for phonons in a suspended crystal in com- 
parison to phonons on a substrate since these phonons 
could travel especially in long nanotubes substantially 
until they diffuse out of the system. This is the reason 
for the experimental finding that large nanotubes show a 
bias response for the G+ mode and the G_ mode in gen- 
eral [2l[. During this travelling the electron-scattered 
optical phonons are then thermalized by scattering with 
acoustic phonons. A true understanding of the different 
behavior of the frequency shift and broadening between 
suspended nanotubes and nanotubes on a substrate re- 
quires a new theory for suspended nanotubes being out 
of the scope of this work. 



VI. SUMMARY 

In this paper we have used the coupled Boltzmann 
equations for electrons and optical phonons to calculate 
the current-voltage characteristic of carbon nanotubes ly- 
ing on a substrate under high bias voltage. First we have 
studied the coupled electron-phonon Boltzmann system. 
By taking into account the electron-phonon relaxation 
time of Ref. [l^ which agrees well with experiments, we 
have determined by numerical fitting the relaxation times 
of the optical phonons in the single- mode relaxation time 
approximation. These are much longer for short nan- 
otubes below 1/zm than for large ones. The result was 
obtained by fitting our numerically determined currcnt- 




FIG. 8: (Color online) The left panel shows the function Fi(x) 
(|38[) . On the right hand side we show the Raman lineshape 
function F2{x) defined in ((39)) . 



voltage curves with the experiments. For short nan- 
otubes, this time did not agree with experimental find- 
ings. In Sect III we went beyond the single time relax- 
ation approximation by taking into account also lower ly- 
ing secondary phonons in the Boltzmann equation which 
leads us to the conclusion that the phonon relaxation 
shows a bottleneck in the sector of acoustic phonons for 
short nanotubes, but not for long nanotubes. We have 
explained this by the fact that due to the phonon veloc- 
ity of the lower-lying acoustic phonons, these phonons 
are redistributed in such a way that at least locally a 
bottleneck is created for short nanotubes. This leads to 
a plug in the relaxation path of the optical phonons. 

In Section IV, we have considered an analytical solu- 
tion of the Boltzmann system where we first linearize the 
electronic equations and use further the assumption of a 
constant phonon distribution function to solve them in 
the electronic sector. We have compared our results for 
the current- voltage characteristic and the phonon distri- 
bution function with the numerical findings of Sect. II, 
Sect. Ill and the experimental results. We find an espe- 
cially good agreement in the high-voltage, small length 
regime. Our analytical theory provides us with the elec- 
tron and phonon distribution functions as a function of 
position and momentum. This opens the possibility to 
calculate the optical phonon broadening and frequency 
shift due to the coupling of the phonon system to the 
bias driven electronic system. For zero bias this coupling 
is the dominant contribution for both quantities 

We have then calculated in Sect. V in a charge-current 
conserved way the broadening and frequency shift of the 
zone-center optical phonons. We find that the phonon- 
level broadening, determined at zero bias voltage mainly 
by the electron-phonon coupling, vanishes at large volt- 
age. The vanishing was explained by a smoothing of 
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the electron distribution function over the Fermi-level at 
non-zero bias voltage. For very large voltages and small 
nanotubcs, we found a positive frequency shift. Note that 
this behavior of a vanishing broadening and positive fre- 
quency shift is in contrast to the self-energy contribution 
of the phonon-phonon interaction to the lower lying hot 
acoustic phonons. This contribution to the self-energy 
led to a negative frequency shift and an additional broad- 
ening at higher bias due to a temperature increase of the 
lower lying acoustic phonons. 

It was in fact just recently found experimentally by 
Raman-scattering, that the level broadening decreases 
and the frequency shift increases for the zone-center op- 
tical phonons 0] for increasing bias voltage. In contrast 
to the very high voltage regime, we found a negative 
frequency shift at moderate bias voltage and large nan- 
otubcs. 



we obtain from ((T|) 



_ flF (t) ^ e _j, 



De 



sinh [ — st I 



with 



De = st cosh | — 



st cosh ( — st ) + 

"or 



sum - — st 

IL 



(A2) 



(A3) 
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where hp{t) is the Fourier-transform of the Fermi- 
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Appendix A: Analytic calculation of the 
current-voltage characteristic 

In the following, we carry out the calculation of the 
current-voltage characteristic and the electron-phonon 
distribution functions for ell ^> Hlj explicitly by using 
(HJ)- dZj) - To solve the system of equations we use the ap- 
proximations ((15)) and (fT^l) . 



fn = h[x -> x — L,x- L x] 



(A4) 



We note that one has to take into account during this 
substitution the absolute value signs in the cosine hyper- 
bolic and sinus hyperbolic arguments in (jA2[) . 



In the following, we evaluate (|A2j) by the help of a 
residuum integration. With 



h = ft° + ff° 



(A5) 



we obtain for the electron distribution function /£ _0 due 
to the residuum at t = 



1. Current- voltage characteristic 

In the following calculation we go in Fourier-space 

fL/n{k L /R,x) = l/{2ir) 2 f dkdtf L/R (t,k)cxp{i(-kx T 
tvphk^m)] and use further the func- 
tion f L/R defined by fL/R(k L / R ,x) = 

l/(2ir) f dtf L/R (t,x)cxp[Ti(tv F hk L/R )}. With the 
help of 



st = \ 1 - 



2ncos(<jjt) + exp(-iojt) 



-I 2 



2n+ 1 



(Al) 



ft=0 

Jl 



[1 - n F {e - eE(x - L))} 



[np(e — eEx) — np(e — eE(x — L))] 



\x-L\ 1 

F 1 4- A 
l sc 1 T 77^ 



(A6) 



f R =0 is given by (jA6|) with the substitution (|A4[) . The 
( ^ singularities are given by the zeros of De. For 
L/lg C ^> 1 we obtain for these zeroes 



st 



st 



for \st\ < 1 , (A7) 
for \st\ > 1 . (A8) 
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With the help of 
A 1 (x 1 ,x 2 ) = t^-t< [l - n F (e - eExx)] 



where H(x) = 7+^(x+l) and 7 is the Euler-Maschcroni 
constant, ^(x) is the digamma function. We show in 
Fig. 6 the function B\{x). 



np (e — eEx\) exp ( — — — |e — eExA 
n uj 



<y~± ln<l + exp 



A 2 {x 1 ,x 2 ) 



2 ! + #■ 



_tt / |e - eExi\ \x 2 \ 

L_\\ UJ 1 ll 



(1 — np (e — eExi)) 



(A9) 



hf (e — eEx±) exp — — — |e — eExA 
n uj 



1 1 



exp 



eEx\ I 



l£2| 



(A10) 



we obtain for /*^° due to the singularities in the regime 
\st\ < 1 JAU) 

= A 1 (x,x-L)+A 1 (x-L,x)+A 2 (x-L,x). (All) 
fpf is given by with the substitution (|A4[) . 



Finally, we have to discuss the contribution of singular- 
ities \st\ > 1 (|A8|) to /*^°. Due to the Fourier-exponent, 
these terms are exponentially suppressed beyond a small 
energy strip of range |Ae|/w > 1 in comparison to terms 
calculated with residua \st\ < 1. From the discussion 
above and the following discussions we obtain that these 
terms are not relevant in the range eU/huj >1, L))!^ 
for the determination of the phonon distribution func- 
tion. We also note without an explicit calculation here 
but which can be shown with similar methods as for the 
current calculation of the \st\ < 1 singularities below that 
the current from the \st\ > 1 singularities is negligible in 
comparison to the current calculated from the singulari- 
ties I si I < 1. This current contribution will be calculated 
explicitly in the following. 

The current / can be calculated by the help of I = 
(4e//i) / de(fn — fp)- This is best done by carrying out 
the integration at the boundaries x = or x = L of the 
nanotube. From (IA6[) and (IA11[) and the correspond- 
ing expressions for /r, wc obtain the following current- 
voltage characteristic 



U 



with 



Bx{x) = 



2 71 



H 



H 



ln(4) 



(A12) 



2H(x) 
(A13) 



2. Phonon distribution function 

Next we calculate the phonon distribution function 
n(x) (fl5|) . We carry out this calculation only in the lead- 
ing order in l^ c /L. In order to simplify our notation wc 
define x = x/L. Next we calculate from (|A6[) . (|Alip 
and the corresponding equation for fp (|A4[) the phonon 
distribution function in the regime ireU 3> huiL/l sc 



a. Phonon distribution function for nell S> fkoL/l BC 

In the following, we only have to determine the k > 
phonons since we find from symmetry arguments that 



n(k, x) = n(—k, —x) 



(A14) 



To simplify our notation we define the following regimes 
for k > 

Tlx : x<^,eU(l-x)> -y->ef/x, (A15) 

1 vphk 
Tl 2 : x > - , eUx > > eU(l - x) , 

1 vphk 1 vphk 
IZ3 : x<—,eUx> or x > —, eU(l — x) > . 

2 2 2 2 

We obtain for the phonon-Boltzmann equation by 
using (|A"6j) . (|ATT|) and ([A4]) 







n(k, x) = Kin(k, x) + K 2 



(A16) 



with 

Tlx : K x 
K 2 : K x 
n 3 : Kx 



-5-1 








±\ 


^sc 




V Top ^ 


1 




-,K 2 = 



- — ) , K 2 

Top J 

s p x(l - x) 

^ph^cp 



sP(l-x) 

^ph^cp 

(A17) 



By taking into account the boundary conditions n(fc, 0) = 
and further v p hT op >C L we obtain for n(fc, x) in Tlx and 

7^3 



72.1 : n ni {k,x) 



s p L 1 



,r, 



(A18) 



Tl 3 : m Z3 (k,x) = — 1 2a(l-a). (A19) 

l!phT e p AxL 



We point out here that according to (|A17[) the expres- 
sions for Kx in (|A18|) and (|A19[1 arc different in the dif- 
ferent regimes (|A15[) . Wc obtain in these regions the 
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same solutions n(k,x) of ([2]) as we would set immedi- 
ately v p h = from the beginning. This is not true for 
n(k,x) in 1Z2- We obtain especially in the regions where 
l sc /L ~ T C p/r p that the neglection of the i> p h term in @ 
is not allowed. Such parameter values lead to the result 
that the phonon distribution function at the boundary of 
the nanotube increases which is in fact seen in Fig. 5. In 
the regime K2 we obtain from (|A16|) . (|A17[) 



and 



K 2 : nn 2 (k,x) = n a (k,x Q (k))e K ^ x o(k)) + 



s p L 



(I + K lX ) - e Kl(x - Xo{k » 



K X L 



VphTep KiL 

(1 + K lXo (k)) 

(A20) 



where x Q (k) = L(l/2 + \l/2-v F hk/2eU\), n (k,x (k)) = 
n-iz 3 (k,xo(k)) for vphk/2 < eU/2 and no{k,xo(k)) = 
hvv F hk/2 > eU/2. 

We now use the abbreviation 



K x = KiL . 



(A21) 



Next, we calculate the phonon density n(x) (fT5"|) . Ac- 
cording to the discussion above, we obtain 

n(x) ft {n Ul (L - x) + n U2 {x))6(x - 1/2) (A22) 
+ (rT Kl (x) + n U2 (L - x))6{l - (x - 1/2)) + 2n n , (x) 

where n-jii are the energy averaged phonon distribu- 
tion functions (|A18|) . (|A19|) and (|A20|) according to (fl"5)) 
where the integration region are restricted to IZi (|A15jl . 
These are given by 



nm (x) 
nu 3 (x) 
and 

nn 2 (x) ft 



1 s p L 1 



2 v ph T ep Ki 
s p L 1 / 1 

VphTep K 1 \2 



^x(l - 2x), 



1 



(A23) 
x(l-x) (A24) 



s p L 1 



1 



VphTep Ki 

2/2 



K\ 



(A25) 

e *i(*-4) 



where we took into account that L/u p hT p ^ 1- Finally, 
we calculate the position averaged phonon distribution 
function n = Jdxn(x)/L. We obtain 



2(n Kl + n n2 + n K J 



with 



n n 3 



s p L 1 1 
VphTep Ki 48 ' 

s p L 1 13 
VphTep K x 192 



(A26) 

(A27) 
(A28) 



niz 2 



s p L 1 

VphTep Ki 



1 

~48 
2 

K* 



1 



(A29) 



By taking into account the definition of (IA21|) . (IA17|) 
we obtain that (|A26|) is the relation which determines the 
energy-position averaged phonon distribution function n 

(EH). 



Let us first assume that Ki 1 in region 1Z2- From 
expression (|A29[) we find that there is no solution in 
this parameter region for ()A26[) since L/v p hT cp ^> 1 and 



V F T C - 



130nm. 



Next we assume that K\ <C — 1 in the region 7?.2 ■ From 
(|A17[) we find then that K\ ft — L/v p hT op in the regions 
TZi. By using (|A26[) we obtain 



_ Us p T op 

n ft - for L ^> 

96 r OD 



Vph Tcp 



(A30) 

For r op — 9.1ps we obtain n = 5.3. This leads to l sc = 
ll.lnm. This value is in excellent agreement with the 
experimentally determined value l sc ft 10 — llnm 0,0]. 
From (|A"22"1) , (|A"23"1) , (lA"24"l) and l|A"2%l> we obtain that the 
phonon distribution function n(x) is zero at the boundary 
of the nanotube. By using the scattering parameters in 
dHJ) we further obtain the validity of (|A30p for L > 
426nm. 

For smaller nanotube lengths one has Ki ~ 0. By 
explicit calculation we obtain for L/v p hT p 3> 1 and K\ ~ 
that n-R^ , n-ji 3 <C n>n 2 . By taking only n-ji 2 into account 
in (|A22[) and (|A26[) we obtain for the energy averaged 
phonon distribution function by Taylor expanding (|A25|) 
with respect to K\ 



n{x) 



^_LQ|i-l/2| 2 -i|£-l/2| 3 ) 



and for (|A26]) with (jA"29jl 



s p L 1 

VphTep 96 



(A32) 



From ([A"3T|) we get an increasing behavior of the phonon 
distribution function at the boundary of the nanotube 
in agreement with Fig. 5. 

Finally, we calculate the phonon distribution function 
n by solving (|A"26l) for n with (|A"27|) . (|A"28f and (|A"29l) 
numerically as a function of the nanotube length L. The 
result is shown in Fig. 7. We obtain that the decrease 
of n for in the direction of small L given by (|A32[) is 
not very large such that we can assume approximately 
the validity of (|A30[) in the whole regime L > 50nm. 
Next, we calculate the phonon distribution function for 
Trell <C HujL/I sc . 
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b. Phonon distribution function for ireU -C huL/l a 



By taking into account ([X6]) . (|A11|) and 
we obtain that n(k,x) 7^ n|j c only in the regime 
\w(vF\k\/2ui)l sc / L\ < 1 where we use that n < 1 as will 
be shown below. In this regime we obtain for k > (|A16|) 
with 



#2 



C0S(7TX) 

1 — cos 2 (7ri) 



1 



(A33) 



cp 



where again (|A14jl holds. We now assume that l sc /L <C 
T e p/r p which is good fulfilled in the large length regime 



where we have (fT4|) T cp /r op 
|7r(» F |*|/2a;)/«/L| < 1 

n(fc, x) s 



0.22 . Then we obtain for 



' O I J 



4 T, 



(A34) 



ep 



We point out that (|A34j) is not valid in the small regime 
where < |e| < eU,nuj. We obtain a position dependent 
numerical prcfactor to (|A34[) in this regime. From (|A34[) 
we are able to calculate n(x) and n given by 



n(x) 



4 T, 



(A35) 



cp 
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